b=fir1(13,0.5)
h=20*log(abs(fft(b)))/log(10)
plot(h)